Army  Research  Laboratory 


Prediction  of 
Electromagnetic  Fields 
Generated  by  Rail  Guns 

by  Ira  Kohiberg 


ARL-CR-148 


September  1 995 


prepared  by 

Kohiberg  Associates,  Inc. 
5375  Duke  Street,  Unit  1 603 
Alexandria,  VA  22304 

under  contract 

DAAL03-91-C-0034 
Subcontract  9341 9 


DTIC 

ELECTE 

OCT  2  4  1995; 

G 


19951023  149 


OnSQUALIWINSPBCrED 


Approved  for  public  release;  distribution  unlimited. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of 
the  Army  position  unless  so  designated  by  other  authorized  documents. 

Citation  of  manufacturer's  or  trade  names  does  not  constitute  an  official 
endorsement  or  approval  of  the  use  thereof. 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the 
originator. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No,  0704-0188 

Pubic  raporling  burdan  for  ttM  coHactiofl  of  information  it  asimatad  to  avaraga  1  hour  par  rasponsa,  including  tha  lima  for  laviawing  inatrucbona,  aaarehing  axiating  data  aourcaa, 
gafoaring  and  maintaining  tha  data  naadad,  and  oomplaling  and  ra^Aawing  tha  eoilaction  of  infornwtion.  Sand  commanta  ragarding  ttaa  buidan  aatimata  or  any  othar  aapact  of  this 
coUaction  of  information,  including  auggaaions  for  radudng  this  burdan,  to  Washington  Haadquartara  Sarvioaa,  Diractorata  for  Infonnation  Opamtiona  and  Raports,  1215  Jaffaraon 
Da^a  Highway,  Suita  1204,  Arfington,  VA  22202-4302,  and  to  tha  Offica  of  Managamant  and  Budgat,  Paparwoik  Raduction  Projact  (0704-01 88),  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  AMMbAnlr;  2.REP0RTDATE  X  REPORT  TYPE  AND  OATES  COVERED 

September  1 995  Final,  from  Jul  93  to  Jan  94 

4.  TITLE  AND  SUBTITLE 

Prediaion  of  Electromagnetic  Fields  Generated  by  Rail  Guns 

S.  FUNDING  NUMBERS 

PE:61102A 

6.  AUraOR(S) 

Ira  Kohiberg 

ARL  POC:  William  O.  Coburn 

7.  PERFORMING  ORGANCATION  NAME(S)  AND  AODRESS<ES} 

Kohiberg  Associates,  Inc. 

5375  Duke  Street,  Unit  1603 

Alexandria,  VA  22304 

•l  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

ARL-CR-148 

«.  SPONSORINCUMONITOflING  AGENCY  NAME(S)  AND  ADDRESS<ES) 

U.S.  Army  Research  Laboratory 

Attn:  AMSRL-WT-ND 

2800  Powder  Mill  Road 

Adelphi,  MD  20783-1197 

10.  SPONSORINGMONTORINQ 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

AMS  code:  P61 1 1 02.00000  Contract  No.:  DAAL03-91  -C-0034 

ARL  PR:  3E25E3  Subcontract  No.:  9341 9 

12a.  DISTRIBUTION/AVAILABIUTY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

12b.  DISTRIBUnON  CODE 

13.  ABSTRACT  (k4aximum  200  words) 

Elearomagnetic  launch  systems  such  as  rail  guns  generate  electromagnetic  signatures  that  contain 
a  spearum  of  frequencies  extending  from  quasi-DC  up  to  tens  of  kilohertz.  However,  electromagnetic 
compatibility  (EMC)  considerations  pertaining  to  the  potential  adverse  effects  of  these  fields  on  nearby 
electronic  equipment  remain  largely  unexplored.  This  report  includes  a  discussion  of  the  theoretical 
models  used  to  predict  the  inductance  gradient,  the  transient  behavior  of  the  currents  produced  in  the 
rail  gun  structure,  the  dynamical  generation  of  the  external  fields,  and  a  comparison  of  the  theoretical 
model  with  experimental  data. 

The  predicted  rail  inductance  gradient  of  =  0.52  pH/m  compares  very  well  with  the  measured 
value  of  0.522  pH/m.  The  existence  of  an  inductance  gradient  efficiency  factor.  Eg,  is  demonstrated, 
with  a  derived  value  of  0.75.  This  produces  an  effective  inductance  gradient  of  L'  =  ^  =  0.39 

pH/m.  Using  this  value  of  L'  leads  to  a  theoretical  exit  muzzle  velocity  of  525  m/s,  which  is  within  5 
percent  of  the  measured  value. 

14.  SUBJECT  TERMS 

Rail  gun,  EM  fields,  magnetic  field 

1Sl  NUMBER  OF  PAGES 

42 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  16.  SECURITY  CLASS! RCATION  19.  SECURITY  CLASSIRCATION 

OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified  Unclassified 

20.  UMTEATION  OF  ABSTRACT 

UL 

NSN  7540-01-280-5500 


Standard  Form  298  (Rav.  2-69} 
PMteribMl  by  ANSI  Sid.  Z39-18 
29S-102 


13.  ABSTRACT  (cont'd) 


The  prediaed  waveshapes  of  the  magnetic  field  components  are  generally  in  good  agreement 
with  observations  for  points  close  to  the  bore  center.  For  radial  distances  greater  than  a  foot,  the 
experiments  show  peak  fields  that  may  exceed  the  predictions  by  a  factor  of  two  to  three.  These  issues 
are  being  investigated. 


AQces(oh"P5r" 

NTIS  CRA&J 
DTIC  TAB 
Unannounced 
Justification 


By _ 

Distribution  / 

Availability  Codes 


Dist 


Avail  and/or 
Special 


Contents 


1.  Introduction  . 5 

2.  Electromagnetic  Model  for  Rail  Gun  . 7 

3.  Inductance  and  Shielding  Considerations . 9 

4.  Physical  Models . 11 

4.1  Force  on  Armature  and  Equations  of  Motion . 11 

4.2  Rail  Gun  Viewed  as  Transmission  Line:  Driving  Inductance  Per  Unit  Length  . i....  13 

5.  Computation  of  Magnetic  Fields  . 17 

6.  Comparison  Between  Theory  and  Experiment  . 21 

7.  Conclusions . 25 

Acknowledgement  . 25 

References  . 26 

Distribution  . 41 

Appendices 

A.  Application  of  Transmission  Line  Theory  to  Rail  Gun . 27 

B.  Equivalence  Between  Electric  and  Magnetic  Fields  in  Two  Dimensions  . 29 

C.  Time  Behavior  of  Magnetic  Field  at  Selected  Field  Points . 33 

Figures 

1.  Rail  gun  and  coordinate  system  . 5 

2.  Rail  gun  electrical  model  . 7 

3.  Rail  gim  in  support  structure . 10 

4.  Model  for  estimating  magnetic  shielding  . . 10 

5.  Top  view  of  rail  gim  . 11 

6.  Computation  of  force  on  armature  . 11 

7.  Transmission  line  model  for  rail  gim . 13 

8.  Magnetic  field  lines  in  upper  right-hand  quadrant . 16 

9.  Inductance  efficiency  as  a  function  of  2a/b . 16 

10.  Geometry  for  computation  of  magnetic  field  contribution  from  armature  . 18 

11.  Geometry  for  computation  of  magnetic  field  contribution  from  rails . 18 

12.  Armature  position  as  a  function  of  time  . 22 

13.  Comparison  between  experimental  and  theoretical  behavior  of  vertical  component  of 

magnetic  field  in  plane  of  symmetry  at  2  in.  above  rails  and  28  in.  from  breech . 22 

14.  Magnetic  field  behavior  at  z  =  28  in.,  y  =  4  in.,  and  x  =  0, 12,  and  24  in . 23 

15.  Armature  and  rail  contributions  to  magnetic  field  intensity,  H(A/m),  at  z  =  30  in. 

from  breech,  at  height  y  =  12  in.  above  rails,  and  in  plane  of  symmetry  {x  =  0) . . 24 

Tables 

1.  EM  model  as  function  of  Dq . 8 

2. fg  as  function  of  b/a  . 15 

3 


1.  Introduction 

In  recent  years  significant  advances  have  been  made  in  the  development  of 
electromagnetic  launchers  (EMLs).  These  devices  can  impart 
hypervelocities  to  macroscopic  size  objects  (e.g.,  greater  than  several 
grams).  Possibilities  even  exist  for  launching  aerospace  platforms  using 
this  technology.  A  rail  gun  is  a  type  of  EML  that  shows  promise  of  being 
implemented  into  operating  systems  in  the  near  future.  Figure  1  shows  a 
generic  model  of  a  rail  gun  that  is  similar  to  the  one  used  in  our 
experiments. 

A  robust  pulse  power  source  (see  sect.  2)  provides  a  burst  of  energy  lasting 
several  milliseconds  that  drives  the  projectile  between  the  rails;  the  projec¬ 
tile  eventually  leaves  the  system  with  a  velocity  in  the  kilometers  per  sec¬ 
ond  range.  The  currents  may  exceed  10^  A,  and  can  generate  very  strong 
magnetic  fields  (H  >  10^  A/ m)  close  to  the  rails,  which  drop  off  with  dis¬ 
tance  from  the  rails.  The  dominant  portion  of  the  electromagnetic  signa¬ 
tures  contains  frequencies  from  dc  to  tens  of  kilohertz.  Depending  on  the 
physical  contact  between  the  projectile  and  rails,  there  may  also  be  inter¬ 
mittent  arcing  at  the  contact  surface.  In  addition,  the  emergence  of  the  pro¬ 
jectile  from  the  rails  is  accompanied  by  arcing,  which  is  alleviated  by  plac¬ 
ing  a  resistive  load  across  the  muzzle. 

As  rail  guns  emerge  from  laboratory  models  toward  implementation  into 
systems,  one  must  consider  the  electromagnetic  compatibility  (EMC)  is¬ 
sues  due  to  the  electromagnetic  fields  generated  by  these  devices  from  the 
enormous  currents  produced  on  the  rails.  The  first  step  in  assessing  EMC 
problems  is  the  characterization  of  these  fields  as  a  function  of  time  and 
distance  from  the  rails.  This  report  describes  the  present  status  of  our  theo¬ 
retical  model  used  to  predict  these  fields,  as  well  as  including  a  limited 
comparison  with  experimental  data.  Recent  papers  by  Coburn  et  al  [1,2] 
address  the  observations  and  experimental  configuration  in  detail. 
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The  incorporation  of  a  rail  gun  into  a  weapon  system  requires  a  quantifica¬ 
tion  of  the  relationship  between  the  power  source,  rail  gun  dynamics,  and 
electromagnetic  fields  generated  by  the  rail  and  armature  currents.  In  or¬ 
der  to  ensure  EMC  with  other  equipment  it  is  necessary  to  characterize  not 
only  the  dominant  low-frequency  (dc  to  tens  of  kilohertz)  fields  generated 
by  the  moving  armature,  but  also  the  higher  frequency  emissions  gener¬ 
ated  by  continuum  arcing  along  the  rails,  and  the  larger  arc  produced  at 
the  muzzle.  For  predicting  the  low-frequency  magnetic  fields,  the  circuit 
approximation  to  the  rails,  which  has  traditionally  been  used  to  study  ar¬ 
mature  acceleration,  appears  to  be  adequate.  On  the  other  hand,  a  modal 
analysis  of  the  rail  system  model  viewed  as  a  transmission  line  is  required 
at  the  higher  frequencies  in  order  to  model  arcing  contributions.  In  addi¬ 
tion,  it  may  be  necessary  to  incorporate  the  effects  of  electromagnetic 
shielding  in  the  evaluation  of  EMC  for  an  actual  system. 

This  report  is  the  first  in  a  series  of  reports  that  eventually  will  provide  the 
ability  to  quantify  EMC  considerations  for  all  frequencies  of  interest.  We 
concentrate  at  present  only  on  the  low-frequency  magnetic  field  generated 
by  the  rail  current  in  a  solid  armature,  for  which  there  is  negligible  arcing. 
It  includes  a  discussion  of  the  pulsed  power  source,  theoretical  consider¬ 
ations  for  predicting  the  fields,  and  a  comparison  between  theory  and 
experiment. 


2.  Electromagnetic  Model  for  Rail  Gun 

Figure  2  shows  a  circuit  model  for  the  rail  gun  that  includes  a  circuit  dia¬ 
gram  for  the  pulsed  power  system.  The  rails  are  connected  to  four  capaci¬ 
tor  banks,  each  composed  of  a  2000-)iF  capacitor  and  a  12-|iH  inductor.  Di¬ 
odes  are  placed  across  the  capacitor  bank  to  prevent  reverse  current.  The 
charging  voltage  is  typically  4  kV.  The  circuit  includes  the  ohmic  resistance 
and  inductance  of  the  source  and  the  connection  to  the  rails.  The  100-0 
load  resistor  serves  as  a  means  to  measure  the  muzzle  voltage. 

The  time  behavior  of  the  current  is  determined  by  combining  the  circuit 
model  of  figure  2  with  the  equations  for  the  motion  of  the  armature.  It  can 
be  shown  that  to  a  good  approximation  the  rail  system  can  be  modelled  as 
a  lossless  transmission  line  with  inductance  and  capacitance  per  unit 
length  defined  as  and  Cr  ,  respectively.  Cr  is  related  to  Lr  by  the 
equation 


and  for  our  application  is  shown  to  be  negligible.  Thus,  capacitive  effects 
are  not  important. 

The  power  supplied  to  the  rails,  VqIo,  is  given  by 

where  the  first  term  on  the  right-hand  side  is  the  magnetic  energy  increase 
per  unit  time,  and  the  second  term,  is  the  power  supplied  to  the  ar¬ 
mature.  Closure  for  the  system  of  equations  is  completed  by  relating  the 
armature  position  Z,<i,  and  the  power  supplied  to  it,  to  the  current.  The  ca¬ 
nonical  equation  for  the  armature  acceleration  is 


m- 


d% 

dt^ 


(3) 


where  m  is  the  mass  of  the  projectile  plus  armature,  and  is  an  effective 
inductance  per  unit  length.  The  selection  of  the  appropriate  value  of  Lg 
has  been  of  longstanding  interest  in  rail  gun  theory.  A  commonly  used 
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Figure  2.  Rail  gun  electrical  model. 


Rails 


value  was  L'^ ;  this,  however,  frequently  led  to  results  that  were  at  vari¬ 
ance  with  experiment.  In  this  report  we  show  that  Lg  can  be  written  as 

L'e=EoL'R,  (4) 

where  Eg  is  an  electromagnetic  efficiency  factor  that  is  always  less  than 
unity. 

When  the  correct  value  of  L  g  is  used,  the  velocity,  ,  is 


0 


Since  the  electrical  power  supplied  to  the  mass,  ViIq,  must  equal  its  rate  of 
increase  in  kinetic  energy,  we  have 

Equations  (2)  to  (6),  when  combined  with  the  circuit  model  for  the  power 
source,  determine  the  current,  voltage,  and  armature  position. 

As  we  show  in  section  4,  Lg  is  determined  from  a  two-step  process.  The 
first  part  is  the  evaluation  of  the  inductance  per  unit  length,  L^,  for  the  rails 
viewed  as  a  parallel  plate  transmission  line.  The  second  part  deals  with  the 
realization  that  only  a  fraction  of  Lr,  which  we  call  the  efficiency  factor,  Eq, 
is  actually  involved  in  the  armature  acceleration  process.  This  leads  to  the 
derivation  of  equation  (4). 

The  electromagnetic  theory  [3-7]  that  leads  to  the  derivation  of  L'^  and  Eq 
also  uniquely  determines  the  electric  and  magnetic  fields,  and  rail  surface 
currents  in  the  x-y  plane  (sect.  4.2).  Strictly  speaking,  the  surface  currents 
are  never  uniformly  distributed  along  the  rail  surfaces,  although  a  uniform 
surface  current  assumption  may  be  acceptable  when  the  ratio  of  rail 
height-to-separation,  alb,  is  large  (see  fig.  6,  p  11). 

For  characteristic  distances  from  the  rail  center,  Dq,  that  are  much  greater 
than  both  a  or  b,  the  detailed  distribution  of  surface  currents  along  the  rails 
is  immaterial  to  the  value  of  the  magnetic  field.  In  this  case  the  approxima¬ 
tion  of  the  rails  as  current  filaments  is  entirely  adequate;  the  filament  ap¬ 
proximation  then  also  applies  to  the  armature.  Table  1  summarizes  the 
choice  of  electromagnetic  field  model  as  a  function  of  Dq  compared  with  a 
and  b.  For  the  uniform  sheet  and  filament  models,  the  magnetic  fields  are 
determined  by  a  straightfoward  application  of  the  Biot-Savart  law.  Since 
experimental  observation  has  been  made  only  in  the  regions  where  Dq  > 
a,  b,  we  discuss  the  analytical  results  for  these  models  in  section  6. 

Table  1.  EM  model  as  EM  model 

function  of  Dq.  Dq  ~b,a  Conformal  mapping 

Do  >  b,  a,  Uniform  sheet  (when  a>b) 

Do  »  b,  a  Filament 
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3.  Inductance  and  Shielding  Considerations 

This  section  addresses  certain  EMC  issues  that  are  important  in  the  inte¬ 
gration  of  rail  guns  into  larger  systems.  Two  basic  requirements  are  in¬ 
volved  here:  (1)  that  the  introduction  of  the  rail  gun  does  not  adversely  af¬ 
fect  the  performance  of  other  systems  and  (2)  that  the  rail  gun  power 
efficiency  be  acceptable.  This  latter  requirement  means  that  for  a  given 
power  source,  a  significant  fraction  of  the  pulse  energy  be  transferred  to 
the  projectile  as  compared  with  increasing  the  magnetic  field  energy 
around  the  rails. 

Figure  3  shows  a  model  of  the  rail  gun  in  the  support  structure  used  in  our 
experiments.  No  special  attempt  was  made  to  shield  the  magnetic  fields.  If, 
however,  a  structure  of  this  type  were  to  be  used  to  shield  nearby  elec¬ 
tronic  equipment,  the  external  frame  would  have  a  relative  permeability, 
fJr,  greater  than  unity.  In  this  case,  not  only  would  the  rail  inductance  gra¬ 
dient  be  changed  from  to  a  new  value,  but  also  the  electromagnetic  ef¬ 
ficiency  would  be  changed.  Theoretical  predictions  for  this  case  have  not 
yet  been  made.  The  essential  challenge  for  EMC  is  to  achieve  the  desired 
balance  between  rail  gun  energy  efficiency  and  shielding. 

On  the  other  hand,  it  is  possible  to  estimate  the  amount  of  shielding  re¬ 
quired  by  approximating  the  rails  as  current  filaments  encased  in  a  cylin¬ 
drical  shell  of  thickness  A  =  r2  -  rj,  as  shown  in  figure  4.  The  solution  for 
the  fields  has  been  worked  out  by  Smythe  [8].  The  magnetic  vector  poten¬ 
tial  has  only  one  component,  A^iz^)  that  is  given  by 


^z-  K 


4n+2 


1-1 


2M-I-1 


«+i 

\2r) 


cos(2  n  -I- 1)0  ,  (7) 


where  fir  is  the  relative  permeability,  and  the  current,  7,  is  a  function  of 
time  determined  from  the  solution  of  the  equations  in  section  2.  The  mag¬ 
netic  field  components  in  cylindrical  coordinates  are  given  by 


^6  = 


(8) 
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Figure  4»  Model  for 
estimating  magnetic 
shielding. 
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4.  Physical  Models 

4.1  Force  on  Armature  and  Equations  of  Motions 

Figure  5  shows  the  top  view  of  the  rail  gun,  while  figure  6  shows  a  blow¬ 
up  of  the  armature  relevant  to  the  calculation  of  the  accelerating  force  in 
the  2-direction.  As  indicated,  figure  6  is  rotated  90°  counterclockwise  from 
figure  5  to  make  the  coordinate  system  compatible  with  the  conformal 
mapping  discussion  of  section  4.2.  We  have  also  indicated  the  existence  of 
surface  currents  in  figure  5  that  are  assumed  to  penetrate  into  the  conduct¬ 
ing  surfaces.  They  do  not,  however,  affect  the  derivation  of  the  armature 
force,  Ffl. 
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Fa  is  determined  by  initially  formulating  the  force  per  unit  volume  in  the 
armature  Jj,  and  then  integrating  over  the  armature  volume.  The  volu¬ 
metric  force  is  given  by 

f^=JxB  ,  (9) 

where  J  is  the  current  density  and  B  =  iiqH  is  the  magnetic  flux  density. 
From  Maxwell's  equations  we  have  the  relationship 

J=VxH.  (10) 

Using  equation  (10)  in  equation  (9)  then  gives  the  total  force. 


mZx+A  _  _ 

(VxH)xH  dxdydz  . 


(11) 


Application  of  the  vector  identity 


(H-V)H  =  iV(HVHx(VxH)  (12) 

converts  the  volume  integral  into  a  surface  integral.  The  z-component  of  F 
is  the  acceleration  force.  Fa  =  Fz,  and  is  given  by 


U. 

S  -b-^ 


h\z  =  Z,)-H\z  =  Z 


+  A)]dxd\/  .  (13) 


The  first  term  in  equation  (13)  is  a  surface  integral  over  all  the  six  surfaces 
of  the  armature  shown  in  figure  6.  For  a  good  conductor,  the  components 
of  H  normal  to  the  armature  surfaces  will  be  very  small,  and  the  first  term 
in  equation  (13)  will  then  be  negligible  compared  to  the  second. 


If  we  assume  that  H^fz  =  Za  +  A)  «  (z  =  Z^),  which  is  consistent  with 
the  generally  assumed  notion  that  there  is  only  a  minimal  magnetic  field  in 
front  of  the  armature,  the  accelerating  force  becomes 


H^z  =  Z^)  dxdy  . 


In  the  next  section  we  show  that 


(14) 


(15) 


where  £o  is  the  previously  defined  efficiency  factor.  The  armature  accelera¬ 
tion  is  then  determined  from  the  equation 


m  IFa^1e±’1^ 

A  ^^2  2  ' 


and  the  power  supplied  to  the  armature  is  F^fdZj^/dt). 


(16) 


4.2  Rail  Gun  Viewed  as  Transmission  Line:  Driving 
Inductance  Per  Unit  Length 

The  purpose  of  this  section  is  to  derive  equation  (15)  using  the  techniques 
of  conformal  mapping.  Accurate  evaluation  of  L'g  is  available  from  C.  E. 
Baum  et  al  [6].  Although  a  rigorous  determination  of  Eq  is  still  in  progress 
(see  app  A),  we  present  an  estimate  of  this  parameter  based  on  an  early  in¬ 
vestigation  presented  by  Assadourian  and  Rimai  [3]. 

Figure  7  shows  the  rail  gun  viewed  as  a  transmission  line.  The  inductance 
per  unit  length,  and  the  rail  capacitance  per  unit  length,  Cr,  are  de¬ 
termined  from  the  solution  of  Maxwell's  equations  for  the  transverse  elec¬ 
tromagnetic  (TEM)  mode.  It  is  this  mode  that  applies  when  transmission 
line  theory  is  used.  The  circuit  model  for  the  rail  gun  is  itself  the  limiting 
form  of  the  transmission  line  version  of  the  rail  gun  system  applied  at  zero 
frequency. 

The  TEM  mode  may  alternatively  be  viewed  as  a  quasi-static  solution 
to  the  class  of  field  problems  solved  by  Laplace's  equation  in  two 
dimensions: 


.d^V 
dx^  dy^ 


=  0  , 


(17) 


where  V  is  the  voltage.  The  capacitance  per  unit  length  for  the  rail  gun  sys¬ 
tem  of  figure  7  is  found  by  first  solving  equation  (17)  for  the  boxmdary  con¬ 
dition  where  the  top  conductor  is  at  potential  Vq  the  bottom  conductor 
at  potential  -Vq-  Surface  charge  density  at  each  conductor,  <Js(x,  y),  is  deter¬ 
mined  from  the  equation 

where  E„  is  the  normal  component  of  the  electric  field  on  the  conducting 
surface,  and  the  total  charge  per  unit  length,  Q,  on  the  top  conductor  is 
found  by  integrating  equation  (18)  over  the  surface.  A  charge  per  unit 
length  of  -Q  resides  on  the  bottom  conductor  in  order  to  close  the  electric 
flux  lines. 


Figure  7. 

Transmission  line 
model  for  rail  gun. 
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L 
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t 
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It  is  important  to  note  that  surface  charge  will  not  only  accrue  on  the  inner 
surfaces  facing  one  another,  but  also  on  the  edges  and  on  the  top  surface. 
In  any  event,  the  capacitance  per  unit  length  is  then  found  from  the 
formula 


r'  —Q- 

^  2Vn 


(19) 


By  expressing  in  the  form  [6] 

^'R  =  ^0^fg  ' 

where  fgisa  geometric  factor  that  depends  only  on  the  ratio,  b/a,  and  then 
using  the  velocity  of  light,  v/. 


gives 


(21) 


(22) 


Table  2  provides  values  of as  a  function  of  (bja)  for  a  width-to-height  ra¬ 
tio  that  may  be  of  interest  in  rail  guns.  For  the  ARL  rail  gun,  (b/a)  =  0.8, 
which  then  gives  =  (0.415)(4  k  x  10“^)  =  0.52  pH/m.  The  method  by 
which  equation  (17)  is  solved  is  that  of  conformal  mapping  [3-7],  which 
has  formed  the  basis  for  analytical  solution  of  Laplace's  equation  in  com¬ 
plex  geometries  such  as  the  rail  gun.  Although  the  field  problem  is  formu¬ 
lated  as  an  electrostatic  problem,  from  which  the  surface  charge  density, 
as(x,  y),  and  the  electrostatic  field,  E  (x,  y),  in  all  x-y  space  is  determined, 
the  properties  of  the  magnetic  field  and  surface  current  distribution  are 
also  readily  available.  It  is  well  known  from  TEM  mode  theory  that  if  a  to¬ 
tal  current  I  is  flowing  on  a  conductor,  the  spatial  distribution  of  current 
density,  js(x,  y),  is 


where  as(x,y)  and  Q  are  computed_for  the  corresponding  electrostatic 
problem.  As  shown  in  appendix  B,  H  (x,  y)  is  given  by 


H  =  ^tJc^xE  ,  (24) 

where  k  ^  is  the  unit  vector  in  the  z-direction. 

The  ratio,  Os(x,  y)/Q,  is  the  normalized  distribution  of  charge  per  unit 
length  that  satisfies  the  path  integration,  and 


cis(x>y)ds  =  'l  , 


(25) 


is  a  function  only  of  b/a,  and  has  dimensions  of  (length)-^  Since  the  elec¬ 
tric  field,_  E,  will  also  be  proportional  to  the  potential  difference,  IVg,  the 
ratio  of  E/Q  will  depend  only  on  the  geometry. 


Table  2.  fg  as  function 
of  b/a  (selected 
results). 


b/a 

4 

b/a 

4 

05 

0.306 

0.9 

0.444 

0.6 

0.346 

1.0 

0.472 

0.7 

0.382 

1.2 

0.522 

0.8 

0.415 

1.4 

0.566 

(Table  source:  Ref  6.) 


In  summary,  the  conformal  mapping  solution  to  Laplace's  equation  ap¬ 
plied  to  the  parallel  plate  rail  gun  uniquely  determines  the  inductance  per 
unit  length,  L'^,  the  distribution  of  current  on  the  rail  surfaces,  and  the  spa¬ 
tial  distribution  of  the  magnetic  field.  Figure  8  shows  a  sample  of  magnetic 
field  lines  in  one  quadrant  of  the  geometry. 

is  defined  by  the  standard  equation 

=  ,  (26) 

where  H  lies  in  the  x-y  plane.  The  integration  of  equation  (26)  is  carried 
out  over  the  entire  x-y  plane,  not  just  the  surface  of  the  armature.  By  com¬ 
bining  equations  (14)  and  (26)  we  derive  the  result 


where  the  efficiency  factor,  Eq,  is  given  by 


+a  r+b 


n 


H^dxiy 


Eo  =  - 


I  J  H^dxdy  +  I  I  H^dxdy 


(28) 


and  r  f  signifies  the  x-y  integration  external  to  the  area  of  the  armature;  this 
is  the  leakage  contribution. 

It  is  possible  to  render  an  estimate  of  Eq  using  an  approximation  devel¬ 
oped  by  Assadourian  and  Rimai  [3].  Their  results,  as  expressed  in  our  ter¬ 
minology,  are  rendered  in  figure  9.  For  the  Army  Research  Laboratory 
(ARL)  rail  gun,  (Talb)  =  2.5,  which  then  gives  an  efficiency  factor  of  about 
75  percent.  The  appropriate  inductance  per  unit  length  to  use  for  accelera¬ 
tion  is  then  L'g  =  EqLr  =  (0.75)  (0.52)  =  0.39  pH/m.  It  will  be  shown  in  sec¬ 
tion  6  that  the  value  of  L'g  provides  good  agreement  with  experimental 
data. 
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Figure  9.  Inductance 


Source:  Reference  3 


5.  Computation  of  Magnetic  Fields 

This  section  addresses  the  computation  of  the  magnetic  fields  for  the  ex¬ 
perimental  configuration  used  in  this  study.  In  this  case  there  is  essentially 
no  shielding  introduced  by  the  support  structure  for  the  rail  gun  system, 
and  contributions  to  the  electromagnetic  fields  arise  from  the  currents 
flowing  on  the  rails  and  armature.  For  each  current  source  the  Biot-Savart 
law  can  be  used  to  compute  the  magnetic  field,  H ,  since  the  wavelengths 
associated  with  the  frequencies  of  interest  are  much  larger  than  the  charac¬ 
teristic  dimensions  of  interest. 

Figures  10  and  11  show  the  geometry  for  computing  the  contributions  to 
the  magnetic  field  intensity  H  from  the  armature  and  rails,  respectively.  If 
we  assume  in  the  general  case  that  the  current  is  distributed  vertically,  so 
that  J(}/)  dy  is  the  differential  current  between  y  and  y  +  dy,  we  have 

(29) 

where  ds  is  a  differential  vector  pointing  in  the  direction  of  the  current,  r 
is  the  distance  from  cE  to  the  observation  point  whose  coordinates  are  a:,  y, 
z,  and  r  is  the  magnitude  of  r.  For  the  armature  ds  =  dx  t,  where  T  is  a 
unit  vector  in  the  Ar-direction,  and  for  the  rails  ds  -  ±dz  k,  where  k  is  a 
unit  vector  in  the  z-direction.  The  "+"  sign  is  used  for  the  rail  whose 
current  is  coming  out  of  the  paper,  and  the  sign  for  inward-directed 
current. 

If  D  =  (X^  +  is  the  transverse  radial  dimension  of  the  field  point,  and 
Do  is  the  maximum  of  either  h  or  w,  there  arise  two  cases  of  practical  inter¬ 
est.  These  are  the  situations  where:  (1)  D  is  slightly  greater  than  Do,  or 
(2)  D  »  Do.  In  the  former  case  the  detailed  structure  of  }(y)  may  be  impor¬ 
tant,  while  for  case  (2)  the  detailed  structure  of  }(y)  is  not  important,  and 
current  filament  approximations  to  the  armature  and  rail  currents  may  be 
used. 

The  distribution  of  surface  current  for  the  rails  is  determined  from  the  con¬ 
formal  mapping  solution  for  the  inductance  per  unit  length  discussed  in 
section  4,  although  for  the  dimensions  of  our  system,  a  uniform  distribu¬ 
tion  might  be  appropriate.  It  is,  however,  not  of  critical  importance  in  this 
report  because  the  observation  points  for  which  the  fields  are  computed 
are  located  in  the  range  D  »  Do,  where  the  filament  approximation  can  be 
used. 

A  rigorous  prediction  of  the  surface  current  distribution  on  the  armature  is 
not  easy  to  obtain.  To  a  first  approximation  one  might  assume  that  its  ver¬ 
tical  distribution  is  the  same  as  that  of  the  rails.  Here  again,  if  D  »  Dq,  the 
current  filament  approximation  for  the  armature  can  also  be  used.  For  fu¬ 
ture  EMC  studies,  where  the  determination  of  close-in  magnetic  fields  may 
be  important  for  implementing  shielding,  it  is  desirable  to  have  available 
an  analytical  expression  for  the  fields  due  to  uniform  sheet  models.  In  this 
way  one  can  obtain  a  first-order  approximation  of  the  importance  of  verti¬ 
cally  distributed  current  sheets. 
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Figure  10,  Geometry 
for  computation  of 
magnetic  field 
contribution  from 
armature. 


y-axis  p(X,  Y) 


Front  view—armature 
coming  out  of  paper 


P(X,  Y) 


Side  view-armature  traveling 
in  positive  Z’direction 


Figure  11.  Geometry 
for  computation  of 
magnetic  field 
contribution  from 
rails. 


(a)  Side  view  of  rails  (b)  Front  view  of  rails 


The  2-component  of  H  due  to  the  armature  is  denoted  as  Hy{zf  the  y- 
component  as  Hj{y^  These  contributions  are  obtained  by  integrating  equa¬ 
tion  (30)  over  the  armature  surface.  Equivalent  expression  for  the  rails  can 
be  obtained  with  appropriate  change  of  coordinates,  but  for  brevity  these 
results  are  not  presented.  For  the  sheet  of  uniform  current  density,  /,  we 
have 


(30a) 
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where 


a‘'’  =  x  +  M;/2  , 
a~  =  x-wll  , 


M,= 


R^-sT^ 


^  rUsT^  ' 


^R--sJa^ 

2  „-.  rr  ' 


R'+Va^ 

R^  =  \/a}  +  2^  +  {Y+hnf  , 

R-  =  ^a^  +  2\{Y-h/2f  , 

N{a)=  I  JP'*'-P1  . 

Va^Z 


(30c) 

(30d) 

(30e) 

(30f) 

(30g) 

(30h) 

(30i) 

(30j) 

(30k) 


|3-  =  arctan|y^ ^4^)  (301) 

Z  =  Z-Z^.  (30m) 

When  X  or  y  are  much  greater  than  w  and  h,  the  foregoing  results  reduce 
to  the  current  filament  approximation,  with  the  total  current  given  by  fh. 
We  then  have 


where 


H- 


1 

(X  +  ic/2)  (X-h;/2)  1 

((X  +  w/lf  +  alf^  ({X  -  w/2f  + 

«J=z^y^ 


(31a) 

(31b) 


(31c) 

(31d) 
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The  magnetic  fields  due  to  the  rails  are  H^x  and  Hj^y,  and  are  given  by 

Hr,  =  -  -  “’/2)  Kk  -  (X + R>n)  KjJ  ,  (32b) 

where 


K 


R  j2| 


b^  =  {X-w/2f  +  Y^  , 


K  ~  — 


I  2  — (—2  —2^/2 

[z  +  b  j  (z  +  b  j 


f  =  (X  +  w/2f  +  Y^  . 


(32c) 

(32d) 

(32e) 

(32f) 
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6.  Comparison  Between  Theory  and  Experiment 

The  primary  emphasis  of  our  study  to  date  has  been  on  the  validation  of 
the  electromagnetic  efficiency  factor  and  the  magnetic  fields  generated  by 
the  rail  gun  systems.  The  acceleration  of  the  armature  is  given  by  equation 
(3).  Using  a  measured  current  that  has  the  waveform 

7o  =  190  X 10^  (exp(-500 1)  -  exp(-1.3  x  10^  f))  (.4)  ,  (33) 

we  can  immediately  integrate  equation  (3)  to  obtain  the  velocity  and  posi¬ 
tion  of  the  armahare  as  a  function  of  time.  Figure  12  shows  a  representative 
comparison  between  the  theoretical  and  experimental  behavior  of  ZaH) 
using  equation  (33)  with  an  11-g  mass,  and  a  theoretically  predicted  value 
of  L'g  =  0.39  |xH/m.  As  observed,  the  agreement  is  good.  The  correspond¬ 
ing  curve  for  velocity,  is  also  in  good  agreement  with  experimental 
results,  yielding  only  a  7-percent  uncertainty  at  the  muzzle. 

Figure  13  shows  a  comparison  between  the  experimental  and  theoretical 
behavior  of  the  vertical  component  of  the  magnetic  field  in  the  plane  of 
symmetry  at  2  in.  above  the  rails  and  28  in.  from  the  breech.  The  theoretical 
result  was  obtained  using  equations  (31),  (32),  and  (33)  and  the  time  behav¬ 
ior  of  ZaH)  determined  from  the  solution  of  equation  (3).  Although  the 
comparison  is  good,  we  are  still  exploring  the  possible  sources  of  error. 

It  is  also  interesting  to  observe  the  space-time  variability  of  the  magnetic 
field  as  a  function  of  position.  A  sample  result  is  shown  in  figure  14,  for 
which  all  field  components  are  indicated  at  the  point  z  =  28  in.,  y  =  4  in., 
and  X  =  0, 12,  and  24  in.,  respectively.  The  observations  are  in  general 
agreement  with  these  results,  although  it  is  found  that  as  we  move  away 
from  the  rails  the  experimental  values  begin  to  exceed  the  theoretical  pre¬ 
dictions  by  more  than  a  factor  of  2.  We  are  exploring  this  matter  further. 

When  the  value  L'  =  Eg  L'^  =  0.39  |xH/m  is  used,  the  position  of  the  arma¬ 
ture,  as  determined  by  the  peaking  characteristic  of  the  magnetic  field  re¬ 
ported  in  the  forward  (z-direction),  is  also  consistent  with  experimental 
observations  [1,  2].  This  is  shown  in  figure  15,  where  theoretical  predic¬ 
tions  of  Hx,  Hy,  and  Hz  using  a  filament  model  for  the  rails  and  armature 
are  rendered  for  z  =  32  in.,  x  =  0,  and  y  =  12  in.  The  predicted  time-to-peak 
of  tp  =  2.4  ms  for  Hz  lies  within  experimental  uncertainty. 

The  experimental  plot  of  Hz  for  the  aforementioned  case  [1]  has  the  same 
waveshape  as  that  shown  in  figure  15  but  is  a  factor  of  2.5  higher.  On  the 
other  hand,  for  observation  points  that  are  close  to  the  rail  there  is  very 
good  agreement  between  the  predicted  and  observed  waveshapes  and 
peak  fields  for  all  magnetic  field  components,  with  uncertainties  ap¬ 
proaching  only  30  percent. 

Appendix  C  provides  a  more  comprehensive  evaluation  of  calculated 
magnetic  fields. 
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7.  Conclusion 


Predictions  for  projectile  acceleration  are  in  good  agreement  with  experi¬ 
mental  observations  for  rail  gun  inductance  and  efficiency  factor.  Muzzle 
exit  velocity  is  found  to  be  within  5  to  7  percent  of  experimental  data.  A 
comparison  between  calculated  and  measured  magnetic  fields  shows  good 
agreement  close  to  the  rails,  with  uncertainties  exceeding  a  factor  of  two  at 
large  distances.  These  sources  of  error  are  being  examined. 

The  large  magnetic  fields  produced  by  these  devices  will  require  shielding 
in  order  to  ensure  compatability  with  other  electronic  equipment.  Shield¬ 
ing  considerations  must  be  considered  as  an  integral  part  of  the  rail  gun 
design  because  it  has  a  direct  effect  on  the  effective  inductance  for  accelera¬ 
tion.  This  in  turn  influences  size  of  the  system  and  terminal  projectile 
velocity. 
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Appendix  A. — Application  of  Transmission  Line  Theory 
to  Rail  Gun 

This  appendix  discusses  the  importance  of  using  a  rigorous  transmission 
line  theory  to  provide  an  accurate  representation  of  the  electromagnetic 
fields  in  the  spatial  region  around  the  rails.  There  are  two  reasons  for  ex¬ 
ecuting  this  calculation.  It  will  provide  an  improved  evaluation  of  the  elec¬ 
tromagnetic  efficiency 


E„  =  - 


I  I  H^dxdy  +  j  I  H^dxdy 


(A-1) 


and  may  also  provide  insight  into  the  coupling  problems  involving  close- 
in  shielding. 

As  we  can  see  from  equation  (A-1)  the  determination  of  Eq  involves  inte¬ 
gration  over  the  leakage  volume  f  f .  The  question  is  how  far  from  the  rail 
bore  center  do  significant  contributions  from  f  f  dx  dy  exist?  An  esti¬ 
mate  of  this  is  available  from  the  early  work  of  Assadourian  and  Rimai  [3], 
which  is  shown  in  figure  A-1.  This  figure  shows  that  virtually  100  percent 
of  the  power  is  contained  over  spatial  dimensions  that  are  of  the  order  of 
the  rail  height.  Although  the  calculations  were  performed  using  a  value  of 
(b/a)  =  0.86  (which  differs  from  the  ARL  value  of  (b/a)  =  0.80),  the  results 
are  not  very  different,  as  shown  in  figure  8. 

Another  important  physical  attribute  of  the  rail  gun  system  is  the  rela¬ 
tively  large  physical  thickness  of  the  rails  as  compared  to  the  separation 
between  them.  Assadourian  and  Rimai  [3]  also  provided  an  estimate  of 
this  effect  which  did  not,  however,  appear  to  introduce  important  modifi¬ 
cations  to  our  conclusions.  In  a  more  refined  calculation  this  aspect  of  the 
problem  should  be  more  critically  addressed.  Figure  A-2  shows  a  com¬ 
parison  in  the  field  lines  between  a  thin  and  thick  strip  above  a  perfectly 
conducting  plane. 

In  addition  to  references  3-7  (in  the  main  body  of  text),  the  author  has  re¬ 
viewed  other  conformal  mapping  approaches  to  the  parallel  plate  trans¬ 
mission  line  system  [9-11].  The  important  issue  is  to  select  the  appropriate 
transformations  in  the  complex  to  simplify  the  computation  of  Eq.  Al¬ 
though  some  work  has  already  been  started  in  this  arena,  the  results  are 
too  premature  to  present  in  this  report. 

In  summary,  a  rigorous  calculation  (e.g.,  conformal  mapping  theory)  is  re¬ 
quired  to  predict  the  fields  in  the  close-in  region  surrounding  the  rails,  and 
to  determine  the  electromagnetic  efficiency,  Eq.  The  detailed  structure  of 
the  currents  on  the  rails  would  appear  to  be  an  important  consideration  in 
the  fields  for  distances  on  the  order  of  several  rail  heights.  At  larger  dis¬ 
tances  from  the  bore  center  the  filament  approximation  should  suffice. 
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(b)  Thick  rail 


Appendix  B.— Equivalence  Between  Electric  and 
Magnetic  Fields  in  Two  Dimensions 

The  theory  of  conformal  mapping  plays  an  important  role  in  determining 
the  spatial  behavior  of  electric  and  magnetic  fields  in  two  dimensions.  Tra¬ 
ditionally,  theoretical  models  are  set  up  as  electrostatic  problems  in  which 
expressions  for  the  electric  fields  follow  directly.  This  appendix  shows 
how  the  electrostatic  results  can  be  easily  used  to  compute  the  spatial  dis¬ 
tribution  of  the  magnetic  fields  when  waves  are  propagated  in  the  trans¬ 
mission  line  mode. 

The  existence  of  the  aforementioned  relationship  has  been  pointed  out  to 
this  author  by  Dr.  Kelvin  Lee,*  although  the  derivation  of  the  result  was 
not  available.  Presumably,  this  result  appears  in  the  open  literature. 
The  equivalence  formula  between  H  and  E  presented  here  is  indepen¬ 
dently  derived,  although  the  result  is  not  new. 

Consider  the  cross  section  for  an  arbitrarily  shaped  conducting  surface 
shown  in  figure  B-1.  Current  and  charge  are  assumed  to  flow  in  the  z- 
direction,  which  in  this  diagram  is  coming  out  of  the  paper.  Let  qsi.  2)  be 
the  charge  per  unit  area  on  the  surface,  and  ;s(  i.,  2)  be  the  z-directed  cur¬ 
rent  per  unit  surfacejength.  The  differential  charge  located  at  position  (. 
on  the  surface  is  qsil,  2)  d£  dz ^nd  the  differential  current  is/s(^,  z)  dt, 
where  dl  is  the  magnitude  oidl. 

Now  let  P(x,  y,  z)  be  any  point  exterior  to  the  conductor  as  shown,  which 
does  not  necessarily  lie  in  the  plane  of  the  conducting  surface  element  d£ 
dz.  The  electric  and  magnetic  fields  at  P  due  to  the  surface  charge  and  cur¬ 
rent  are  respectively  given  by 


(B-1) 

(B-2) 


where  r  is  the  position  vector  from  the  surface  element  to  P,k  is  a  unit 
vector  in  the  z-direction,  and  r  is  the  magnitude  of  P . 

In  the  realm  where  the  concepts  of  capacitance  and  inductance  apply, 
which  is  the  transmission  line  mode  of  propagation,  the  spatial  distribu- 


Private  Communication,  K.  S.  H.  Lee,  Kaman  Sciences,  Corp.,  Dikewood  Division,  Santa  Monica,  CA  (September 
1993). 
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Figure  B-1,  Charge 
and  current  on 
conducting  surface. 


Conducting  surface 


tions  of  surface  charge  and  surface  current  are  similar.  We  can  therefore 
write 


qjj^,z)  =  Q{z)r]{lz)  ,  (B-3) 

j^,z)  =  Kz)r\(^ ,z)  ,  (B-4) 

where  Q(z)  is  the  charge  per  unit  length  at  position  z,  J(z)  is  the  current  at  z, 
and  rjlP,z}  is  a  normalized  spatial  distribution  function  that  satisfies  the 
condition 


J  r|(?,  z)d^  =  I  .  (B-5) 

The  integration  of  equation  (B-5)  is  completely  around  the  conducting  sur¬ 
face  at  fixed  z. 

When  the  conditions  for  conformal  mapping  apply,  the  spatial  variations 
with  z  are  neglected.  Thus,  Q,  I,  and  77  are  no  longer  functions  of  z.  In  this 
domain  we  form  the  vector  cross  product,  k  x  E,  using  equation  (B-1). 
The  result  is 


Using  equation  (B-4)  in  equation  (B-2)  we  deduce 


(B-6) 


(B-7) 
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By  comparing  equations  (B-6)  and  (B-7)  and  using  the  formula  B  =  PqH  , 
we  derive  the  desired  result 


H  =  -Q£Qk><E  .  (B-8) 

Equation  (B-8)  is  used  in  conjunction  with  conformal  mapping  applied  to 
electrostatics  in  the  following  way.  Suppose  the  electrostatic  problem  is  set 
in  accordance  with  the  discussion  of  section  4.2.  The  charge,  Q,  appearing 
on  either  rail,  as  well  as  the  electric  field,  E,  are  both  proportional  to  the 
initially  assigned  arbitrary  voltage,  Vq.  The  ratio,  E  /Q,  is  an  intrinsic 
property  of  the  system,  that  is,  independent  of  Vq-  By  simply  specifying  a 
total  current,  I,  flowing  on  either  rail,  we  can  readily  determine  the  mag¬ 
netic  intensity,  H,  from  equation  (B-8). 


Appendix  C . — Time  Behavior  of  Magnetic  Field  at 
Selected  Field  Points 


This  appendix  presents  results  for  the  magnetic  field  behavior  at  selected 
points  to  show  the  variability  with  distance  along  the  rails  (fig.  C-1  to  C-6), 
variability  with  height  above  the  rails  in  the  plane  of  symmetry  (fig.  C-7  to 
C-1 1 ),  and  the  variability  with  displacement  off  the  plane  of  symmetry  (fig. 
C-8  and  fig.  C-1 2  to  C-1 5).  The  legend  in  these  figures  is  explained  below. 


In  these  figures  the  following  coordinates  are  used: 


Rails 


Used  in  text 

Used  in  figures 

fiAz 

AY 

AZ 

HRy 

Hr, 

RY 

RX 

AY  +  RY 

Figure  C-1.  Magnetic 
field  behavior  at 
x:  =  0  in,,  y  =  2  in., 
and  r  =  3  in. 
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Figure  C-4.  Magnetic 
field  behavior  at:c  =  0, 
y  =  4  in.,  and  z  =  32  in. 


Figure  C-5,  Magnetic 
field  behavior  at  a:  = 
0,  y  =  22  in.,  and  z  = 
32  in. 


Magnetic  field  intensity  (x  1 0®)  (A/m)  Magnetic  field  intensity  (x  1 0®)  (/Vm) 


Figure  C-8.  Magnetic 
field  behavior  !Ltx  = 
0,  y  =  4  in,,  and  z  =  28 
in. 


Figure  C-9.  Magnetic 
field  behavior  at  x  = 
0,  y  =  8  in,,  and  z  =  28 
in. 


Magnetic  field  intensity  (x  10^)  (A/m)  Magnetic  field  intensity  (x  10^)  (A/m) 
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